1 



FRAGMENTATION IN KINEMATICALLY COLD DISKS 

D. Huber and D. Pfenniger 

Geneva Observatory, CH-1290 Sauverny, Switzerland 



o 
o 

(N 

c 

a 



> 
(N 



Oh 
i 

O 



13 



Abstract 

Gravity is scale free. Thus gravity may form similar 
structures in self-gravitating systems on different scales. 
Indeed, observations of the interstellar medium, spiral disks 
and cosmic structures, reveal similar characteristics. The 
structures in these systems are very lumpy and inhomo- 
geneous. Moreover some of these structures do not seem 
to be of random nature, but obey certain power laws. 

Models of slightly dissipative self-gravitating disks show 
how such inhomogcneous structures can be maintained on 
the kpc-scale. The basic physical processes in these mod- 
els are self-gravity, dissipation and differential rotation. 
In order to explore the structures resulting from these 
processes, local simulations of self-gravitating disks are 
performed in 2D and 3D. We observe persistent patterns, 
formed by transient structures, whose intensity and mor- 
phological characteristic depend on the dissipation rate. 
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1. Introduction 

Molecular clouds reveal hierarchical structures, obeying 
power laws over several orders of magnitude in scale. Ob- 
servations suggest, that the hierarchical structure of kine- 
matically cold media is not only present in Milky Way 
molecular clouds, but is also found in other systems and 



on larger scales. Vogelaar & Wakker (1994) found, e.g. 
perimeter-area correlations in high- velocity clouds. Power- 
law power spectra were found in the Small and the Large 
Magellanic Cloud by |Stanimirovic ct al (1999) and Elme 



green et al. (2001), respectively. Furthermore, measure- 



ments of the HI distribution in galaxies of the M81 cluster 



reveal fractal structures on the galaxy disc scale (Westp- 
fahl ct al. 1999). The matter on cosmic scales is also hi 



erarchically organized. A common feature of the ISM and 
the cosmic structure is, that the matter distribution can 
be characterized by a comparable fractal dimension. All 
this suggests, that a general scale free factor is mainly 
responsible for the matter distribution and the dynamics 
of cosmic structures, disks and molecular clouds. There is 
only one factor being able to have a dominant influence on 
all these scales, namely gravity. The local shearing sheet 



experiments of Toomre & Kalnajs (1991, hereafter TK) 
show that gravity in combination with shearing and dissi- 
pation can develop long-range correlations and maintain 
the system out of equilibrium. In order to study in detail 
the structures resulting from these processes on the kpc 
scale, we perform local simulations of self-gravitating disks 
in 2D and 3D. The third dimension becomes important as 
soon as a strong matter clumping causes a tight coupling 
of self-gravitating forces in the 3D equations of motion. 

2. Local Model 

2.1. Principle 

In local models of disks, everything inside a box of a given 
size is simulated and more distant regions in the plane are 
represented by replicas of the local box. In such a model 
the orbital motion of the particles is determined by Hill's 



approximation of Newton's equation of motion (Hill 1878) 
In 3D they read 



x - 2Q y = 4O A a; + F x (x,y,z) 
y + 2tt x = F y (x,y,z) 



(1) 



= —vz 



+ F z (x,y,z) , 



where A = — ^TZo(dQ / <1TZ)tz is the Oort constant of dif- 
ferential rotation and v is the vertical epicycle frequency. 
F x , F y and F z are local forces due to self-gravitating par- 
ticles. 

The local system is periodic in the y-direction and iso- 
lated in the z-direction. If we use an afhne coordinate sys- 
tem whose pitch angle changes with time then the sys- 
tem is also in the y-direction periodic (Hubcr & Pfcnnigei 
19991 , 2001). Thus we can calculate the forces with the 



convolution method using the FFT algorithm. Thereby 
the computation time is reduced to be proportional to 
N c log(-/V c ), where iV c is the number of cells, taken here as 
proportional to the number of particles. 

2.2. Weak Friction 

To counteract the particle dynamical heating, TK pro- 
posed to add an ad-hoc friction term playing the role of 
the dissipative factors at work in the interstellar gas. Fol- 
lowing TK we include linear friction terms as well. The 
friction terms should be weak in order to keep a quasi- 
Hamiltonian system. Indeed at the kpc scale the physics 
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Figure 1. The evolution of the particle positions seen from above the galaxy plane. The structures result from a 3D 




is still dominated by gravitational dynamics and its con- 
crete behavior should be weakly dependent on the partic- 
ular dissipative factors. 

The linear friction terms — C x x and —C z z added to 
the radial resp. vertical forces (F x , F z ) control the particle 
motions via Equation [D. There is no azimuthal friction in 



order to be consistent with a global angular momentum 
conservation. 

2.3. Initial Conditions 

We are interested in the secular time behavior of the galaxy 
disk. Thus we perform simulations for t = 10 galactic ro- 



3 




tations. In the initial state at t = the particles are dis- 
tributed uniformly in the x-y-plane. In the ^-direction the 
particle distribution follows an isothermal law 

p oc sech 2 (z/z ) , (2) 

where p is the density and zq is the disk scale height. The 
velocities at t = are determined by the shear 

* = , y = ~2A x , z = (3) 

and the Schwarzschild velocity ellipsoid 
3.36GE Q 



<T X K 



(4) 



y 2n 

where the Safronov-Toomre stability criterion is Q > 1. 
3. Results 

We present here some preliminary results of an extended 



study to this subject flHuber fc Pfenniger 2001 



3.1. 3D SHEARING BOXES 



The structures resulting from shearing box simulations 
depend a lot on the relative strength of the competing 
gravitational and dissipation processes. Gravitational in- 
stabilities lead to a conversion of directed kinetic energy 
(shear-flow) into random thermal motion. In this way the 
disk is heated up. If the dissipation is too weak then ini- 
tially arised structures can not be maintained and smear 

out quickly. If the dissipation is increased, a filamentary N(R) oc M{R) oc R D ^ 



structure can be maintained in a statistical equilibrium. If 
we continue to increase the dissipation the filaments be- 
come denser and denser and clumps in filaments may be 
formed. If finally the dissipation dominates completely the 
heating process hot clumps, collecting almost all the mat- 
ter of the simulation zone are formed out of the filaments. 
The change of the structure morphology for an increas- 
ing dissipation is showed in Figures 1-3. For convenience 
we call the dissipation strengths, leading to the presented 
structures, "weak", "middle" and "strong". The relative 
cooling times corresponding to these dissipation strengths 
are: r cool4 : r coo i !2 : r cooli3 : r osc = 16 : 15 : 14 : 1, where 
r osc is the period of the unforced epicyclic motion. 

The simulations are carried out with 131040 parti- 
cles, corresponding to a surface number density of n = 
3640A^ rit . The dynamical range is 1.8 dex. Only each forth 
particle is shown. 

3.2. Mass-Size Relation 

In order to characterize the structures resulting from the 
shearing box experiments, we determine the mass-size re- 
lation. We choose a representative set of particles and 
count for each particle the number of neighboring particles 
N(R) inside a certain radius. If we repeat this for other 
values of R we can find the structure dimension D(R) via 



(5) 



d\n(R) 

where R denotes the scale. The mass-size relation is then 

(6) 
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If the structure dimension is independent of the scale, D = 
Df, i.e., if D is constant or oscillates around a mean value, 
then the mass-size relation is a power law, 

M oc R D f . (7) 

If furthermore Df is non-integer, the structure is fractal. If 
however the structure dimension depends on the scale D = 
D(R), the structure dimension may simply be regarded 
as a statistical measure describing the dumpiness at the 
corresponding scale. 

In order to asses the general relevance of the under- 
lying physical processes we check the structure for self- 
similarity, i.e., we check if D = Df for a certain scale 
range. However, one has to take into account that the 
structures result form a finite simulation, modeling a fi- 
nite physical system. Thus the system can not be fractal 
beyond an upper and lower cutoff. An upper cutoff is given 
by the scale at which the system become periodic. A lower 
scale limit is due to the finite resolution of the simulation 
mesh. 

Figure |4| shows the structure dimension of the mass dis- 
tribution presented in Figures 1-3. The longer term evolu- 
tion of the structure may be superimposed by fluctuations 
on time-scales of the order of ~ 1/2 T rot , where T ro t is the 
time for a rotation around the galaxy center. In order to 
eliminate these fluctuations we indicate in this paper mean 
values of the structure dimension D determined during the 
last two rotations. 

The structure dimensions are not constant over the 
corresponding dynamical range and are thus not fractal. 
However the structure dimension resulting from the sim- 
ulation with the "middle" dissipation has a structure di- 
mension 1.5 < D < 1.8 over the whole dynamical range 
and remains smaller than 2 also on scales where the disk 
thickness becomes important. Thus the corresponding mat- 
ter distribution can approximately be described by a power- 
law for the considered scale range. The structure dimen- 
sion has a minimum at R = 0.25 A cr j t . The increase of 
the structure dimension at R > 0.25 A cr it | may then be 
due to the lower and upper cutoff. A larger dynamical 
range may thus flatten the curve depicting the structure 
dimension. This supposition is supported by low resolu- 
tion simulations showing a steeper increase of D beyond 
the minimum. 
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4. Conclusions 

The structure resulting from the local simulations of self- 
gravitating disks can be homogeneous, filamentary or clumpy 
depending on the relative strength of the competing gravi- 
tational and dissipation processes. As long as the structure 
is mainly filamentary self-gravitation and dissipation en- 
sure a statistical equilibrium, i.e., persistent patterns con- 
sisting of transient structures are formed. If the dissipative 
processes begin to dominate the evolution, the structures 
turn from filamentary to clumpy. During the subsequent 



Figure 4- The structure dimension D of the mass distri- 
butions showed in Figures 1-3 as a function of the scale 
R. The corresponding simulations were carried out with 
a "weak", a "middle" and a "strong" dissipation, respec- 
tively. The solid line depicts the initial mass distribution. 
On large scales this state represents a 2D matter distribu- 
tion, whereas on small scales it tends to D = 3. 



evolution the clumps become hotter and more massive. In 
general clumpy structures do not evolve towards a statisti- 
cal equilibrium. However 2D simulations with a dynamical 
range of 2.5 dex show that it is also possible to establish 
persistent patterns formed by clumps. 

A larger dynamical range produce in general a flat- 
ter structure dimension curve. However the scale range of 
the simulations is still too small to draw final conclusions 
about self-similarity in open, self-gravitating systems. 
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